Propagation of local decohering action in distributed quantum systems 
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We study propagation of the decohering influence caused 
by a local measurement performed on a distributed quantum 
system. As an example, the gas of bosons forming a Bose- 
Einstein condensate is considered. We demonstrate that the 
local decohering perturbation exerted on the measured region 
propagates over the system in the form of a decoherence wave, 
whose dynamics is governed by elementary excitations of the 
system. We argue that the post-measurement evolution of the 
system (determined by elementary excitations) is of impor- 
tance for transfer of decoherence, while the initial collapse of 
the wave function has negligible impact on the regions which 
are not directly affected by the measurement. 
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INTRODUCTION 

The theory of quantum measurement begun in the 
1920s still remains an active topic of interest (see, e.g. 
Ref. Q and references therein). According to von Neu- 
mann's theory of measurement unitary evolution of 
a system prepared initially in a pure quantum state is 
interrupted by an instant decohering action of the mea- 
suring apparatus, so that the density matrix describing 
an ensemble of such systems changes radically (it ceases 
to be a projection operator) and entropy rises. This view 
has been shown to describe rather accurately the con- 
sequences of an act of measurement, but the dynamics 
of the measurement process itself is lacking. The con- 
temporary theory of quantum measurements, which pro- 
vides much deeper analysis of the measurement process, 
is based on the concept of decoherence ||. To be mea- 
sured, the system has to interact with its environment, 
which consists of a large number of degrees of freedom. 
The Hilbert space of the system becomes divided into 
subspaces corresponding to the same eigenvalue of the 
system-environment interaction Hamiltonian. As a result 
of this interaction, coherence between different subspaces 
is quickly lost, and after the measurement the system 
appears in a mixed state. The concept of decoherence 
turned out to be successful in many areas of fundamen- 
tal physics, such as the study of macroscopic quantum 
effects [|J and consistent histories interpretation of quan- 
tum mechanics g), so that investigation of this process 
and related effects is of considerable importance. 

At present, decoherence and its consequences for point- 
like quantum systems have been studied in detail (for 
review, see Ref. ||), but distributed quantum systems 
have received significantly less attention. Mostly, linear 



systems have been investigated, where separation into 
noninteracting modes is possible, and each mode is con- 
sidered as an independent oscillator However, this 
approach is difficult to apply to sufficiently nonlinear sys- 
tems (e.g., spin systems, or the Bose-Einstein condensate 
as described by the Gross-Pitaevskii equation) possess- 
ing localized soliton-likc excitations. For systems where 
localized excitations prevail, dealing explicitely with real- 
space coordinates could be a more suitable strategy. 

A real-space description of decoherence in distributed 
systems is a very general and complicated issue. In 
this paper, we consider only one aspect of the problem, 
namely, how local properties of different regions in a dis- 
tributed quantum system are affected by a local measure- 
ment, that acts only on some part of the system. Indeed, 
different regions in the system are not isolated from each 
other, and correlations between them exist (or can build 
up). Therefore, in spite of the fact that a local measure- 
ment initially affects only one region, other regions can 
"acquire knowledge" that some part of the system has 
been measured. In this paper we explicitely show that 
the decohering influence of the local measurement prop- 
agates through the system in the form of a decoherence 
wave. Dynamics of the decoherence wave is governed by 
elementary excitations, while the effect of entanglement 
is very small for macroscopically large systems. 

The consideration presented here can be applied to 
other similar situations, so that a decoherence wave prop- 
agating with a characteristic velocity of excitations is 
likely to be quite common. This phenomenon, being a 
notable part of any real measurement, is of fundamental 
interest. Moreover, propagation of decoherence can be 
also of importance for the design of quantum computers. 
Such a computer is a system of interacting quantum en- 
tities, representing quantum bits (qubits). Fault-tolerant 
quantum computations involve measurements performed 
on some qubits and it is important to know how such 
measurement may affect other qubits B. Moreover, de- 
coherence is introduced by a dissipative environment of 
qubits, so that analysis of decoherence propagation may 
lead to strategies to minimize influences detrimental to 
performance of the computer. 

In this paper we consider a Bose-Einstein condensate 
of an ideal or weakly non-ideal gas of bosons, which con- 
stitutes a good example of a distributed system in a pure 
quantum state. It can be implemented in reality as a 
gas of trapped atoms cooled down to very low tempera- 
tures H . Suppose we measure the number of particles in 
some region of space. If two such measurements are done 
simultaneously at two different parts of the trap we ob- 
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tain the trivial result corresponding to the ground-state 
wavefunction of the condensate. But if the second mea- 
surement is carried out after some delay then the result is 
different and provides information about the propagation 
of the perturbation induced by the first measurement. 

The situation considered here is related to the prob- 
lem of broken gauge symmetry and existence of a rela- 
tive phase of two interfering condensates JhJ] , which has 
been extensively discussed recently. If we have a con- 
densate with a definite number of particles, its phase 
is spread uniformly between and 2tt, while a definite 
phase requires a non-conservation of the number of par- 
ticles in the condensate. It has been shown that a well- 
defined phase (evidenced experimentally by appearance 
of the interference fringes) builds up in the course of the 
measurement (atoms detection) , due to increasing uncer- 
tainty in the number of particles in each of the interfering 
condensates: each detected atom may well belong to ci- 
ther of them. For the circumstances considered in this 
paper, we have a similar situation: the local phase of the 
condensate is the same in every region. Identity of the 
phase throughout the condensate is due to uncertainty 
in the local number of the particles inside each region. 
However, when the number of particles in some region is 
determined by a local measurement, the phase coherence 
in the condensate as a whole is partially destroyed, what 
leads to observable consequences, propagation of the de- 
coherence wave in the system. Note that decoherence 
wave is the same both for the condensate with definite 
number of particles (with uncertain global phase, the case 
of non-interacting bosons) and for the condensate with 
definite global phase (but with uncertain number of con- 
densed particles, the case of weakly interacting bosons): 
the results for the latter case transform exactly to the 
results for the former as interaction goes to zero. 

We describe the dynamics of the condensate in a linear 
approximation, i.e. we use the approximation of nonin- 
teracting quasiparticles to study a weakly non-ideal Bose- 
gas. In so doing, we loose the ability to investigate some 
interesting nonlinear effects, but we gain in clarity of pre- 
sentation: it is reasonable to start from a simplified (and 
not totally unrealistic) case to emphasize the main idea. 

We do not specify the way of measuring the local den- 
sity of condensate, and the dynamics of the measurement 
process is not considered here. Analysis of a specific ex- 
perimental scheme is a distinct problem, requiring sepa- 
rate study, while here we focus on the post-measurement 
evolution of the condensate. In principle, the local den- 
sity of the Bose-condensate can be measured by placing 
some probe into the trap, which interacts with the con- 
densate so that an entangled state is formed 

\X) = Y,C n \n)®\a n ) (1) 

where \n) is the state of condensate with the number of 
particles n in the measured region, and \a n ) is the state 
of the probe. If the probe interacts with a large num- 
ber of environmental degrees of freedom, so that \a n ) 



are the eigenstates corresponding to different eigenvalues 
of the probe-environment interaction Hamiltonian, then 
the coherence between different probe states is being lost, 
and the condensate's state also becomes an incoherent 
mixture of different states \n). If the probe (and, cor- 
respondingly, the condensate) decoheres quickly enough 
(as is usually the case) we can consider the measurement 
as instantaneous and safely use von Neumann's theory 
to describe the condensate's state immediately after the 
measurement. 

Although the situation considered above is in many 
respects too idealized to apply rigorously to a real ex- 
periment, it is detailed enough to capture the essential 
processes of concern in this paper. 

PROPAGATION OF DECOHERENCE IN 
BOSE-EINSTEIN CONDENSATE 

To study quantitatively the effect of decoherence prop- 
agation, let us consider first an ideal Bose-gas confined 
by external fields and described by the Hamiltonian 

H = Y,E ll ala„ (2) 

A" 

where aj^ and are the boson creation and annihila- 
tion operators. are the one-particle energies, and we 
denote the corresponding one-particle wavefunctions as 
<Pn{r), where n = stands for the ground state having 
minimal energy E$ = 0. Then, the ground-state eigen- 
function of the system of M bosons can be written as 




where |0) is the vacuum state. For simplicity, we can 
consider the trap as being divided into a large number 
N c of small cells each having the volume Vo (it can be 
considered as the volume directly affected by the measur- 
ing apparatus), satisfying the relation Vo -C V, where V 
is the total volume of the trap. Then, the coordinate r is 
understood as a discrete quantity (the number of a cell). 
This is similar to a general practice in solid-state theory, 
where Vo is analogous to the volume of an elementary cell 
of the crystal jy]]. Note that in so doing, the number of 
one-particle states taken into account becomes equal to 
N c , which is finite, though very large. This corresponds 
to the fact that the number of states inside the first Bril- 
louin zone equals to the number of lattice cells. 

At the instant t = we perform measurement of the 
number of bosons in the cell r = 0. This observable is 
represented by the operator N — cv(0)a(0), where 

a ( r ) = '52 < Pn( T ) a »- ( 4 ) 

is the boson field operator. Eigenvalues of the operator N 
are n = 0, 1, 2... and, suppose, the measurement has given 
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us one of them. According to von Neumann's theory, 
it corresponds to the action of the operator W n on the 
system, where 



It; 

/ — exp [i(p(n — N)] 



(5) 



is a projector onto the state with the number of particles 
n in the measured region. The operator W n has the value 
equal to unity on this state and it has zero value on all 
others states. Further development of the system is to 
be described by the density matrix of the system U(t), 
since the measurement interrupts unitary evolution and 
casts the system into mixed quantum state. According 
to the standard theory of measurement , the density 
matrix at the time t is 



U(t) = exp {-iHt)W n U in Wl exp (iHt), 



(6) 



where U- ln = \^>)(^\ is the density matrix before the mea- 
surement. 

To trace propagation of decohcrcncc in the system, we 
study evolution of the one-particle density matrix 



p(r,r',i) = Tr[[/(i)at(r>(r)] 



(7) 



This quantity describes local properties of the Bosc- 
Einstein condensate; in particular, the average number of 
particles resulting from the second measurement, which 
is performed at the point r at the instant t, is given by 
the value p(r, r, t). 

To simplify calculations, we use the fact that the total 
number of particles is large, M ^> 1, so that operators 
cto and a\ acting on the state |\&) can be replaced by 
the number \f~M with relative accuracy 1 / \[~M\ this is 
a standard ap pro ximation in the theory of Bose-Einstein 
condensation |14|| . Therefore, Eq. (^) can be rewritten 
as 



a(r) = y/n B (r) +a(r), a(r) = ^ ^(r)a^ (8) 

where ns(r) = MiPq(t) is the average number of con- 
densate particles contained in the volume Vq at the cell 
r. The expression for the one-particle density matrix can 
be written as 



p(r,r',i) = ^p„(r,r',i) 



(9) 



p„(r,r',t) = (*\Wlai(r',t)a(r,t)W n \V), 



where a(r,i) = exp (iHt)a(r) exp (— iHt). The operator 
product in Eq. (^) is to be ordered normally, i.e. it is to 
be rewritten in such a way that all stand to the left of 
all a in each term of the Taylor series expansion. In so 
doing, we take into account that 



[a(r,t),at( )] =^^(r)^(0)e- 



iE u t 



g(r,t). (10) 



Note that for a system containing a large number of par- 
ticles M » 1, the function g(r, t) can be replaced by the 
Green's function 



G(r,t) = 2^('K(0K 



iE u t 



(11) 



with accuracy of order of 1/JVf, since G(r,t) = g(v,t) 
</?o(r)</3g(0). Performing the calculations, we obtain 



(r,r',<) 



Pn 



vWfr) - G(r, t)y/n^ 



(12) 



yJn B {r')-G*(r',t)^ 
+p„_m G(r, t)G*(r',t), 



where uq — ub(0), and p n = e~ n ° n^/fal) is the Poisson 
distribution function. Summation over n can be per- 
formed explicitly, yielding 



P (r, r', t) = y/n B (T)n B (T') - G*(r', t)^/n B (r)n (13) 
-G(r, t)^/n B {r')n + 2n G*(r', t)G(r, t). 

This result shows that the measurement made at the 
point r =0 produces a decohering perturbation which 
propagates over the trap in the form of a decoherence 
wave, and this propagation is governed by the Green's 
function G(r,t). It can be explicitely demonstrated by 
considering an example of the gas consisting of free Bose- 
particles of mass m. The corresponding Green's function 
at the distances r ^> Vq^ 3 and times t 3> mV^ 3 /h is |T^ 1 



G(r,t) = K)( 



m \ 3 / 2 



/ imr 2 \ 



2mht) GXP V 2?r^i 



(14) 



Local density of the condensate after the measurement is 
given by the value 



(15) 



-2n B V 



-) 



3/2 



2nht) 



cos 



2nht 



where n B = M/V is density of the condensate before the 
measurement, which is independent on position r. This is 
an observable effect, which, in principle, can be detected 
experimentally. 

The entropy of the system, being initially zero, after 
the measurement is 

oo 

S = - Tr [U(t) In U(t)} = - Pn ^Pn > 0, (16) 

n=0 

which is a clear indication of the decohering effect of mea- 
surement. The increase of entropy of condensate as a 
whole happens only at the instant of measurement and 
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further evolution, being unitary, keeps it constant (de- 
coherence only propagates in the system from one re- 
gion to another) . Note that local entropy (in contrast to 
the one-particle density matrix, where the decoherence 
propagation is clearly seen) can be hardly used to track 
the decoherence wave. The value of the local entropy 
is nonzero even in the initial pure state, while the to- 
tal entropy of the system is zero. It happens because of 
"negative entropy" stored in the form of correlations be- 
tween different parts of the condensate (for more detailed 
discussion see Ref. |l3|]). 

The results obtained can be qualitatively interpreted 
as follows. The measurement performed at r = leads to 
localization of some number of particles within the cell 
r = 0. The localized particles acquire rather large mo- 

1 /3 

menta, of order h/V ; the average number of such par- 
ticles is no = n B (0). Immediately after being localized, 
these particles start to propagate over the trap, and their 
propagation is governed by the Green's function ([u]) . Be- 
cause of indistinguishability of particles in the trap, we 
can not say that these are "the same" particles which 
were measured at r — 0, so that the effect we consider is 
not a physical motion of some separate particles in the 
trap, but is the propagation of the decohering influence 
of the measurement through the system. 

An interesting feature of the decoherence propagation 
can be illustrated by the gas of bosons trapped in a 
parabolic external potential, so that each particle is rep- 
resented by an isotropic harmonic oscillator of eigenfre- 

1 /3 

quency fi. In this case, provided that r 3> V and 
V <C (h/£l) 3 / 2 ~ V, the Green's function has the form 
E3 



G(r,t) = V 



2iriTi sin Qt 



3/2 fiVlY 2 

exp — -cot fit] (17) 
\ 27m ' 



where the particles are assumed to have unitary mass. 
This function is periodic in time with the period 27r/fi. 
Therefore, the decoherence propagation is also periodic 
in time with the same period. In the general case of Bosc- 
gas trapped in a finite volume, the decoherence propaga- 
tion becomes a quasiperiodic process, according to Eq. 
©■ 

And, last but not least, decoherence propagation is a 
wave process, possessing both amplitude and phase. Ex- 
istence of coherent waves in the system without quantum 
coherence is not unusual, the same property is shared, 
e.g., by the sound wave propagating in the classical fluid. 
Therefore, in principle, an interference of two decoher- 
ence waves is possible. 

Above, we have considered the system of noninteract- 
ing bosons. Now, let us investigate the case of weakly 
interacting particles, i.e. a weakly non-ideal Bose-gas 
contained in a trap of large volume V. We assume no 
external potential acting on the particles, so that the 
one-particle states are simple plane waves 



exp (ikr) , 



(18) 



where the normalization reflects the fact that the trap 
is divided into cells of volume Vo *C V. This system is 
described by the Hamiltonian 



(19) 



1 

2V 



ki+ks^k'j+k^ 



k 'i) a ik' a k'. a k 2 ak! 



where u(k) is the Fourier transform of the interaction po- 
tential (which is assumed to be repulsive). Since the in- 
teraction is small, new Bose operators can be introduced 
according to Bogoliubov transformation 



«k = £k coshxk + £L k sinhxk, 



(20) 



"_ k - CkSinhxk + Cl k coshxk, 
with the parameters Xk defined as 

v(k)n B 



tanh 2%k 



-Ek + v(k)n B ' 



(21) 



where n B is the average number of particles belonging 
to Bose-Einstein condensate contained in the volume Vq. 
Provided that the interaction is small (or the gas den- 
sity M/V is small), almost all particles belong to the 
condensate, so we can take n B = MVq/V with relative 
accuracy of order of y/v{Q)M/V Q. By using the Bo- 
goliubov transformation, we pass to the ideal gas of new 
excitations with the dispersion law 



Wk = JE 2 _ + 2E k v(k)n B . 



(22) 



Again, we consider dynamical behavior of the one- 
particle density matrix. The calculation procedure re- 
mains essentially the same as for the ideal Bose-gas. In 
so doing, we obtain the result: 



Pn(r,r',t) 



n B 



& 2 



(n!) 2 dz n dz' n 
x[l + (*'-l)G*(r',f)] 

x exp [nsAT(z, z')] > 

where the following notations were used, 

X(z, z') = B{zz' - 1) + (1 - B)(z + z' 
+A[(z-l) 2 + (z>-l) 2 ], 



{[l + (z-l)G(r,t)} (23) 



-2) 



(24) 
(25) 



.4 



1 ^ V ( k ) n B 



2V 



2V ^ 
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and G(r, t) is the Green's function of the weakly inter- 
acting Bose-gas: 



G(r,t) =]Texp(ikr) 



(26) 



, , .E k + v(k)n B . 
x < cosojki — i sinwkt 



Again, we see that the decoherence wave propagating in 
the system follows the dynamics of the Green's function 
(pff). Dynamic behavior of G(r,t) at large times t and 
large distances r can be analyzed by the method of sta- 
tionary phase jL5j . According to this method, the value 
of the function G(r, t) at the point r at the instant t is de- 
termined by those excitations which have a group veloc- 
ity u(k) = duj^/dk obeying the requirement u(k) = r/t. 
The excitations with large wavevectors k are subject to 
considerable damping jl4j , so that at large distances only 
the undamped long- wavelength excitations determine the 
dynamics of the Green's function. These excitations rep- 
resent sound propagating in the Bose-gas with the veloc- 
ity c = y^nBv(0)/m, so the decoherence wave in a system 
of weakly interacting bosons propagates with the sound 
velocity c. 

This result can be interpreted in the same way as the 
decoherence wave in an ideal Bose-gas. The measure- 
ment affects the particles situated at r = 0. Due to the 
interparticle interaction, the decohering perturbation is 
transferred to other regions of the system. The decoher- 
ence transfer is governed by the undamped excitations 
present in the system, i.e. by the long-wavelength exci- 
tations traveling with the sound velocity c. 



DISCUSSION 

Summarizing, we have studied the decohering influence 
of a local measurement performed on a distributed quan- 
tum system. We show that the decohering perturbation 
exerted on the measured region propagates over the sys- 
tem by forming a decoherence wave, whose dynamics is 
determined by the Green's function of the system. This 
result, although not totally unexpected, is not as trivial 
as it might seem, since decoherence is a rather peculiar 
effect, and the decohering impact of a measurement can 
be quite different from other physical influences (see, e.g. 
the discussion in Ref. ]T^|). 

The usual scenario for few-particle systems is based 
on the Einstein-Podolsky-Rosen (EPR) situation [[l7| of 
strong entanglement, when, e.g. two particles with spins 
1/2 form a singlet state 



the transfer of decoherence is instant: the second spin im- 
mediately occurs in the state 1 1). In distributed systems 
this effect is also present: the wave function of the sys- 
tem collapses immediately after the measurement. But 
the impact of the collapse upon the one-particle density 
matrix (and even s-particle density matrix, for s <C M) 
is practically unobservable for the system of macroscopic 
size (where M 3> 1): the change in the density matrix 
element p(r,r',t) immediately after the measurement is 
of order of n /M (provided, of course, that r,r' ^ 0), 
and the same is true for the fc-particle density matrix 
p(ri, . . . r,t; r[, . . . r' k ) if k <C M . This result is rather ob- 
vious: localization of the number hq of particles in some 
cell can not affect noticeably other cells if the total num- 
ber of particles is macroscopically large. Therefore, the 
post-measurement evolution of the system, which is gov- 
erned by the Green's function, becomes important since 
it provides much more noticeable changes in the density 
matrix elements: in Eq. (15) the term corresponding to 
the decoherence wave does not go to zero as M — ► oo. 
Obviously, it happens because in the EPR-like situation 
the entanglement is very "stiff" , so that each state of one 
particle determines completely the state of the other. But 
in the many-particle system there is no one-to-one corre- 
spondence, since the total number of degrees of freedom is 
much larger than the number of degrees of freedom fixed 
during the measurement. This difference is the reason 
for the different dynamics of decoherence propagation. 

Finally, we remark that another aspect of decoherence 
in distributed systems has been studied within the con- 
text of decoherent quantum histories |jl8|,[l9) . Although 
the effects studied there, as well as systems considered 
and methods used, are different from those investigated 
here, it is interesting to note that local properties of dis- 
tributed quantum systems are often "intrinsically" deco- 
herent |l9| if a coarse enough description is used. For the 
effects considered here, sufficient coarse graining leads 
to averaging of the oscillating Green's function over the 
spatial scale of several oscillations, so that the details 
of the decoherence wave becomes negligible. Therefore, 
the intrinsic structure of the decoherence wave can be 
distinguished only at fine scales, where coherence of the 
Green's function holds. 

This work was partially carried out at the Ames Lab- 
oratory, which is operated for the U. S. Department of 
Energy by Iowa State University under Contract No. W- 
7405-82 and was supported by the Director for Energy 
Research, Office of Basic Energy Sciences of the U. S. 
Department of Energy. 



= ^(IU)-UT)). 



(27) 



If the first spin has been measured, and as a result of this 
measurement has been cast in the state 1 f) (here again we 
use von Neumann's theory of instant measurement), then 
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